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Abstract 



Accurate determinations of the physical scale of a lattice action are required 
to check scaling and take the continuum limit. We present a high statistics 
C*") . study of the static potential for the SU(3) Wilson gauge action on coarse lattices 

(5.54 < (3 < 6.0). Using an improved analysis procedure we determine the string 
tension and the Sommer scale (and related quantities) to 1% accuracy, including 
all systematic errors. Combining our results with earlier ones on finer lattices, 
we present parameterizations of these quantities that should be accurate to about 
1% for 5.6 < (3 < 6.5. We estimate the A-parameter of quenched QCD to be 
A m = 247(16) MeV. 



1 Introduction 



One of the technical but very important parts of a lattice simulation is the determination 
of the physical scale of the lattice. Only with accurate scale determinations can one 
check if the ratio of an observable to the physical scale has the theoretically expected 
scaling violations (e.g. order a 2 for the Wilson gauge action). This is required for the 
accurate and reliable continuum extrapolation of observables, and for their conversion 
into physical units to compare with (or predict) experiment. 

Accurate scale determinations have been given new impetus by the recent progress 
in the Symanzik improvement program (see e.g. [ffl, |2|). In particular, the rather non- 
trivial, non-perturbative determination || (cf. also ||) of the 0(a) coefficient in the 
Sheikholeslami-Wohlert (SW) quark action coupled to the Wilson gauge action has to be 
verified by accurate scaling checks. Only then would we feel confident to proceed with 
other steps in the improvement program. 

It is known that the rho mass is a sensitive indicator of scaling violations, in particular 
of the 0(a) violations of chiral symmetry that are generically present in a Wilson- type 
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quark action like the SW action. The rho mass is an ideal quantity to use in scaling 
checks. Therefore we would rather not use it to set the scale, as is sometimes done. 
Instead we have to determine some other quantity sufficiently accurately to set the scale. 

In pure gauge theory (or the quenched approximation of full QCD) the string tension 
a obtained from the static potential is such a quantity — in principle. In practice 
technical problems have often interfered with accurate determinations of this quantity. 
This is especially true on coarse lattices, j3 < 6.0 for the Wilson gauge action^] where, 
otherwise, scaling checks could be performed most accurately. 

To our surprise we found that no accurate string tensions are known for the Wilson 
gauge action for j3 < 6.0. The newest results we are aware of were presented in refs. [|], f| . 
The errors are relatively large, and combined with string tensions on finer lattices the 
results do not really seem to lie on a smooth curve. 

We have therefore decided to perform high statistics simulations of the static potential 
for the Wilson gauge action with (3 < 6.0, to determine the string tension and Sommer- 
type scales [|| (the latter are also well-defined for full QCD). In sect. 2 we describe 
procedures to accurately calculate these quantities from the static potential. Sect. 3 
gives the details of our simulations and the numerical results. In sect. 4 we present an 
accurate parameterization of our results for all f3 in the range 5.6 < /? < 6.5. In sect. 5 we 
discuss the issue of scaling and extracting the A-parameter of QCD from the measured 
string tensions. We conclude in sect. 6. 



2 Scale Determination with the Static Potential 



A standard method of fitting lattice Monte Carlo data for the static potential between 
an infinitely heavy quark and anti-quark is to use an ansatz of the form || 



V(r) 



Vq + or 



+ / 



(2.1) 



Here [£] denotes the tree-level lattice Coulomb term, viz, 



4lT 



d 3 k 

(2tt) 3 



cos(k-r) Z) o(0,k) , 



(2.2) 



where Doo(k) is the time-time component of the gluon propagator for the gauge action 
in question. Note that I is an effective parameter introduced to model lattice artifacts 
beyond tree-level. 

Clearly, eq. ( |2.lD is not a fundamental ansatz, since neither the running of the QCD 
coupling nor the lattice artifacts are taken into account properly (for a more ambitious 
attempt see [|l0|). This is a problem if one tries to extract Aqcd or other (intermediate- 
to) short-distance aspects of the running coupling from potential data via an ansatz of 
the above or similar form. 



1 We use the standard notation /3 = 6/g 2 to parameterize the bare gauge coupling of the SU(3) Wilson 
plaquette action. 
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We are here not interested in short-distance properties of the potential. Instead, we 
would like to extract a physical scale from the (intermediate- to) long-distance part of 
the potential. 

Let us first consider determining the string tension a, which, according to the flux 
tube picture and MC studies, is a well-defined quantity in pure gauge theory. When 
trying to determine only the string tension, the problems of the ansatz (|2.1| ) at short 
distances become largely irrelevant. In practice we will always leave out the first few 
lattice points and check if the string tension stabilizes as we include in the fit only points 
with r > r m j n , with increasing r m j n . 

According to the string picture the potential in pure gauge theory behaves for large 
distances as 

V(r) = V + ar - — + .... (2.3) 

r 

If the transverse fluctuations of the QCD flux tube are described by the simplest bosonic 
string theory, then the "IR charge" takes on the universal value erR = 7r/12 in four space- 
time dimensions [pj| . Like others, we have observed that even with very good statistics 
(and smearing, etc) the intermediate- and long-distance part of lattice potential data 
can indeed always be described by an ansatz of the form ( |2.lD with e = tt/12. 

To determine the string tensions we therefore perform three types of fits: 



4-parameter fits to (2.1) 



3-parameter fits to (2.1) with e = tt/12 fixed 



2-parameter fits to (2.1) with e = tt/12 and I = fixed. 



As r m ; n increases for the various fits, they should reach consistent values for the string 
tension. This asymptotic value is expected to emerge first for the 4-parameter fit, next 
for the 3- and finally for the 2-parameter fit. The demand that all three fits give a 
consistent value should allow an accurate and reliable determination of the string tension. 
We should remark that similar ideas were used in the string tension determinations by 



the Bielefeld group [12] 



Although the method just sketched works quite well, as we will see, there are clearly 
drawbacks to using the string tension to set the scale. First of all, of course, in full 
QCD the string tension is not well-defined, due to string breaking. Secondly, due to 
its long-distance nature, one needs the potential in a region where errors are becoming 
large. 

To avoid these problems Sommer || suggested to define an intermediate distance 
scale r c via the force between a heavy quark and anti-quark, 



V\r c ) = r 2 c F(r c ) = c (2.4) 



for some suitably chosen real number c. The question is how to implement this idea in 
practice to extract r c from lattice data. 



3 



One approach is to define the force and its error via discrete (correlated) differences 
of pairs of points, 

= WM-vm\ 

|ri - r 2 | 

where 77 is chosen || to eliminate tree- level lattice artifacts, 

i _ nil -lill (26) 



rf |ri - r 2 | 

In practice one finds that the force values from different pairs of points around r c still do 
not fall on a smooth curve, due to remaining lattice artifacts. When interpolating the 
force to obtain r c one therefore ends up with relatively large errors. 



One can do better (see e.g. [ 13 [ ) by defining rj via the corrected (more continuum-like) 
potential data, 



V c (r) = V MC (r) + (e-i) 

as 



1 



1 



(2.7) 



fVz) = . + 4 = ^- V ^\ (2i 
or rj \t\ — r 2 | 

where a, e and I are taken from a fit to the ansatz ( |2.1[) . By (slight) abuse of notation we 
are also using V c (r) = Vq + or — e/r on the left side of eq. (|2.8|) to denote the continuum 



part of (2.1). 



Taking the last idea to its logical conclusion we will use the following method to 



determine r c : Perform a correlated fit of the potential data to an ansatz of the form (2.1) 
- but only locally, using data between some r m \ n and r max close to r c . We then obtain 
r c from the continuum part of the ansatz as 



a 



(2.9) 



Since we are not performing a global fit of the potential data the precise from of the 
(continuum part of the) ansatz should not be crucial.0 We can estimate the systematic 
error of r c by varying r m i n and r max , and by fixing e — n/\1. Similarly, we can get an 
idea of the effect of lattice artifacts beyond tree-level by comparing free I with I = fits. 

The nice feature of this strategy is that we can avoid taking the numerical derivatives 
required for the force, while still using only local information about the potential and 
being able to take lattice artifacts beyond tree-level into account. The fits required are of 
the same kind that one performs to determine the string tension; one just has to consider 
smaller intervals [r m ; n ,r max ] around r c . 



2 One would therefore expect the general 4-parameter fit to be somewhat over-parameterized. Indeed, 
wc find that for r c fits e and a are strongly (anti-) correlated. 
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p 


Volume 


L (fm) 


^onfgs 


6.0 


16 3 -32 


1.5 


4000 


5.85 


16 3 -32 


2.0 


8000 


5.7 


16 3 -32 


2.6 


4000 


5.6 


12 4 


2.6 


2000 




8 4 


1.7 


4630 


5.54 


12 4 


2.9 


2000 



Table 1: Volume, spatial extent and number of configurations used in our simulations. 



3 Simulations and Results 

We performed simulations at five couplings f3 < 6.0 for the SU(3) Wilson plaquette 
action. Our updating algorithm alternates microcanonial over-relaxation and heatbath 
steps acting on SU(2) subgroups fl4j (typically in a 3:1 ratio per sweep). In all cases 
we generated several thousand (essentially) independent configurations on lattices whose 
spatial extent is at least 1.5 fm; cf. table |]. We should remark that a subset of the 32-16 3 
configurations at (3 = 5.7 and 5.85 are also being used for spectrum calculations with the 
non-perturbatively improved SW quark action ||. 

The static potential is obtained in the standard manner from (smeared) Wilson loops 
W(r, t) via the effective local masses 

V t (r) = In W ^\. (3.1) 

To reach plateaus for small t we perform APE smearing E^ ] on the spatial link variables. 
This amounts to adding the spatial staples to a suitable multiple of the given link, 

U k (x) aU k (x) + U l( x ) U *( x + U -i( x + * + (3.2) 

and projecting back onto SU(3). This procedure is then iterated n times. 

Note that (modulo details of the SU(3) projection) the above procedure is essentially 
equivalent to a more intuitive "covariant gaussian smearing" in terms of the transverse 
laplacian Aj_, 

U k (x) -> (l + eA ± ) n U k (x) « e enA ^U k (x). (3.3) 

In momentum space the exponential becomes a gaussian, at least in the free case. The 
relation between e and a is0 e ~ 1/(4 + a). 

This way of writing the smearing is useful, because it shows that as long as e is 
small the actual amount of smearing will only depend on the product en. Dimensional 

3 We do not write an exact equality here, because of possibly different implementations of the 
projection, and the question of whether one uses the initial or the "updated" gauge field in the covariant 
laplacian. 
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analysis then also shows that the same amount of smearing, in physical units, corresponds 
to scaling en oc 1/a 2 as one varies the lattice spacing. One therefore has to optimize 
the smearing only at one coupling and can then use this scaling to find a (near-) optimal 
smearing at other lattice spacings.^] We find that en ~ 4 is a good smearing for (3 = 6.0. 
Even for large r ground state overlaps are then at least 85%. 

It is important to have off-axis r to accurately determine the string tension and r c . 
We measured the potential at multiples of the lattice vectors (1,0,0), (1,1,0), (1,1,1) 
and (2,1,0)|] 

Our fitting procedure was sketched in sect. 2. For a given fit, i.e. for a given ansatz, 
t in eq. (jO|), as well as r m i n and r max , we use a single elimination jackknife procedure 
(after suitable binning of the data) to obtain the errors of all quantities of interest. The 
final value and error we quote is obtained as follows. For each fit ansatz (with 2, 3, or 4 
parameters, cf. sect. 2) and t we consider the best few fits (five or less, say), that satisfy 
a minimal "goodness" criterion. Our "goodness" function is typically the product of the 
confidence level and the number of degrees of freedom of the fit. We then combine all 
the "good" fits with a weight that depends on the goodness and the error of the quantity 
in question. 

This procedure is semi-automated to eliminate subjective bias. It is hard to 
completely automate this process. For example, it can happen (e.g. if the smearing is not 
optimal) that fits from neighboring pairs of time slices seem roughly equally "good" , but 
differ significantly in their estimate of the quantity under consideration. In such cases it 
seems preferable to invoke human intelligence and experience to decide what weight to 
assign to the different fits. 

Overall our analysis procedure is conservative, because we do not simply take the 
one "best" fit according to some necessarily somewhat subjective criterion. By always 
taking several fits into account we find that our final numbers and error bars are quite 
stable under changes of the goodness function or of the cuts on the minimal required 
goodness of a fit. Our final error is never smaller than that of any single good fit. In 
case one fit is much "better" than the others, our error and central value will be close 
to the values of this fit; in general our error bar is significantly larger, and our central 
value somewhere between that of the different fits. 

Using this procedure we have measured the string tension and several r c . The 
"standard" choice of r c corresponds to c = 1.65. Slightly inconsistently, but in 
accord with accepted conventions, we will denote this choice by vq. It corresponds 
to tq ~ 0.50 fm (perhaps a bit less); which is around the region where phenomenological 
potentials are best determined. On fine lattices it also corresponds to the region where 
the potential can be very accurately measured. On lattices coarser than 0.1 fm this 
choice becomes increasingly impractical, however, because ro is then quite small in lattice 
units. We have therefore decided to measure and r$ in addition to ro- We note that 
r$ ~ 1.0 fm. 

4 The situation is more complicated for improved gauge actions, due to the contributions from "ghost 
states" to correlation functions at small times. 

5 For the simulations on 8 4 lattices at (3 — 5.6 we measured several more r values. This case is also 
exceptional in that we used an updating algorithm that alternates over-relaxation and Metropolis steps. 
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0.086 



0.086 



fS = 5.85 
t = 4/5 



b 



X 4-par fit 
□ 3-par fit 

o 3-par fit 



p = 5.85 
t = 5/6 



X 4-par fit 
□ 3-par fit 
o 8-par fit 



Figure 1: String tension for (5 = 5.85 as a function of r min for different "good" fits and 
using time slices t = 4/5 and t = 5/6. The solid and dashed lines denote our final value 
and its error, respectively. 



p 




r /o 


r 4 /a 


r 6 /a 


5.54 


0.5727(52) 


2.054(13) 


3.381(21) 


4.187(27) 


5.6 


0.5064(28) 


2.344(8) 


3.814(15) 


4.718(23) 


5.7 


0.3879(39) 


2.990(24) 


4.926(35) 


6.128(42) 


5.85 


0.2874(7) 


4.103(12) 


6.733(17) 


8.340(24) 


6.0 


0.2189(9) 


5.369(9) 


8.831(21) 


10.89(3) 



Table 2: Our results for the string tension and Sommer-type scales. 



A final remark about the r c determinations. Recall that in this case we vary r m [ n 
and r max around r c . We find, perhaps not too surprisingly, that the dependence on r max 
is always small, in most cases negligible. 

Some of our results are illustrated in figures [l] - [3|. Note in particular the stability 
of our determination of ro in figure ||, even for values of r m - m larger than tq itself! Such 
fits are not included in our final analysis, but we regard this as a sign of the stable and 
conservative nature of our analysis strategy. Our final results are collected in table ||. The 
relatively large errors in the (3 = 5.7 case are presumably due to non-optimal smearing; 
we did not follow the guideline described earlier (yet). 

As explained above, we believe that the values quoted include all systematic errors 
- except possibly for one source we have not discussed so far: finite volume effects. 
Previous studies for /3 > 6.0 indicate that a spatial extent of about 1.5 fm is sufficient 
to suppress finite volume errors in the string tension to the level of about one percent or 
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/? = 6.00 



5.35 



X 4-par fit, t=3/4 
+ 4-par fit, t=4/5 
H 4-par fit, t=5/6 



X 3-par fit, t=3/4 
o 3-par fit, t=4/5 
□ 3-par fit, t=5/6 



5.30 



Figure 2: Sommer scale tq for (3 = 6.0 as a function of r m j n for various "good" fits. The 
solid and dashed lines denote our final value and its error, respectively. 



3.1 



3.0 



= 5.70 



cd 
o 



2.9 



3.8 



¥ ^ 



o 3-par fit, t-3/4 
3-par fit, t=4/5 
□ 3-par fit, t=5/6 
x 4-par fit, t=4/5 



Figure 3: Sommer scale vq for (3 = 5.70 as a function of r m i n for various "good" fits. 
The solid and dashed lines denote our final value and its error, respectively. Note that 
the 3-parameter fits with t = 3/4 have r m \ a > vq. These fits were not included in our 
final error estimate; they are shown only for illustration. 
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Scale 


8 4 (1.7 fm) 


12 4 (2.6 fm) 


Change 


a^/a 


0.5004(23) 


0.5064(28) 


+1.2% ± 0.7% 


r /a 


2.355(7) 


2.344(8) 


-0.5% ± 0.5% 


r 4 /a 


3.867(10) 


3.814(15) 


-1.4% ±0.5% 


r 6 /a 


4.787(14) 


4.718(23) 


-1.5% ±0.6% 



Table 3: Check of finite volume effects for /3 = 5.6. 



less [16]. To explicitly check this for our determinations, we have repeated the simulations 
at /3 = 5.6 on 8 4 in addition to 12 lattices. The comparison is shown in table [|. It would 
seem that within our errors there are significant finite volume effects on the 1 — 1.5% 
level, except for ro , which seems to be much less affected (not too surprisingly) . The sign 
of the change of the string tension is as expected (from earlier studies) for finite volume 
effects. That the r c move in the opposite direction is also expected. 

In any case, for coarse lattices there is no problem in using a sufficiently large volume 
to suppress finite volume effects to significantly less than one percent. This is true for 
the values in table ||, except perhaps for /3 = 6.0. We will return to this point in the next 
section. 



4 Parameterizing the Scale 

The string tension is one of the quantities that has traditionally been used to check 
perturbative scaling. It is known that, on a quantitative level, the string tension follows 
neither two- or three-loop scaling, at least for the lattice spacings we are interested in (we 
will see this explicitly below). Allton suggested that this is due to a 2 ,g 2 a 2 , . . . , a 4 , . . . 
lattice artifacts, that are unavoidable if one determines the string tension with the Wilson 
gauge action. 

He therefore proposed jl7|] to fit the string tension to an ansatz of the form 
(a^)(g) = f{g 2 ){l + c 2 a(g) 2 + c A a{gf + ...)/±, a(g) = f Jyl 6)) > t 4 - 1 ) 



in terms of the fit parameters A/y/a, c 2 , C4, Note that a{g) has been normalized to 

1 for /3 = 6.0; we will also use this convention in other coupling schemes considered later. 
In the above f(g 2 ) is the universal two-loop scaling function of SU(3) gauge theory, 

/(^MWPW^), 6o = JiL 6l = _HL (4.2, 

Recall that if scaling holds with a scaling function / in a scheme S, then As = qf{g%{<i)) 
is a constant, the A-parameter of the scheme. 

On theoretical grounds one might want to add g 2 a(g) 2 , . . . corrections to the terms 
shown in eq. (f4.1|), but in practice this does not turn out to be useful (however, we 
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p 


V 


ro/a 


Volume 


Reference 


5.54 


0.5727(52) 


2.054(13) 


12 4 




This work 


5.6 


0.5064(28) 


2.344(8) 


12 4 




This work 


5.7 


0.3879(39) 


2.990(24) 


16 3 


32 


This work 


5.85 


0.2874 7) 


4.103(12) 


16 3 


32 


This work 


6.0 


0.2189(9) 


5.369(9) 


16 3 


32 


This work 




0.2182 16) 




16 4 




Perantonis/Michael |l9| 




0.2184(19) 


5.35(+2)(-3) 


16 4 




SESAM || 




0.2209 23 




32 4 




Bali/Schilling/Hoeber [ 


m 




0.2154(50) 


5.47(11) 


16 3 


48 


UKQCD El 


6.2 


0.1610(9) 




32 4 




Bali/Schilling/Hoeber [ 


20 






0.1604(11) 


7.37(3) 


32 4 




SESAM El 




0.1608(23) 


7.29(17) 


24 3 


48 


UKQCD ||| 


6.4 


0.1214(12) 


9.89(16) 


32 4 




Bali/Schilling g§ 




0.1218(28) 


9.75(17) 


32 3 


64 


UKQCD H 


6.5 


0.1068(9) 


11.23(21) 


36 4 




UKQCD ||, H 



Table 4: Selected world results for the string tension and tq. 



have checked that the A-parameter obtained from fits is stable when using such terms 
in Q). 

In table ^ we have collected some results for the string tension and tq. We are not 



showing data for /3 = 6.6 and 6.8, since they might be affected by finite volume errors IS]. 
Looking at the different volumes employed for the data shown in table f|, it seems quite 
clear that finite volume errors in the string tension are on the one percent or less level 
for all (3 (cf. also our remarks in sect. 3 and the discussion in the references quoted in 
the table). 

These data are well fitted by an ansatz of the form (fOj) using the two-loop scaling 
function for the bare coupling |] We obtain for the 12 degree of freedom fit of the string 
tension with all data in the table, 

A , x 

-= = 0.01338(12) , 



a 

c 2 = 0.195(16) , 

c 4 = 0.0562(45), (4.3) 

with a confidence level of Q = 0.19. Adding an d(<?) 6 term in the ansatz gives a Q of 
0.24. For the three-parameter fit of ro we obtain Q = 0.12. 

When leaving out the j3 = 5.54 point from the fits, the confidence levels jump up 
to Q = 0.9 for three-parameter fits of both ^fa and ro. This might indicate that the 

6 In the next section we will discuss fits in terms of other couplings, that we consider to be more 
physical. Here our aim is the purely pragmatic one of providing an explicit, analytic parameterization 
of the data, for which, it turns out, the bare coupling is just as suited as any other one. 
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determination of these quantities at (5 = 5.54 is increasingly tainted by lattice artifacts 
that are not modeled with sufficient accuracy by our ansatz. We find some signs of this 
in our potential fits. 

We would like to present a parameterization of the data that reproduces the central 
values within one percent and one standard deviation, whichever is smaller (considering 
the tiny error bars in ^/a at /3 = 5.85 we allow a difference of, say, 1.5 standard deviations 
there). This is difficult to achieve with the /3 = 5.54 point, so we will do so only for the 
interval 5.6 < (3 < 6.5. The following parameterizations satisfy these constraints: 

(ay/o)(g) = f{g 2 ) (1 + 0.2731 a{g) 2 - 0.01545 d(#) 4 + 0.01975 a{gf )/0.01364 , 
(a/r )(g) = f{g 2 ) (1 + 0.2W6a(g) 2 + 0.05492 a(g) A )/0.01596 . (4.4) 

The representation of the string tension is shown in figure together with the data from 
table H and results obtained earlier for (3 < 6.0. 



5 Scaling and A^jg 

One can consider taking the above much more seriously than just as a way of 
parameterizing the data. That is, one can try to extract the A-parameter of QCD from 
it. The statistical uncertainties of the fit parameters are very small, about 1% for A/i/a 
and slightly larger for Ar^. However, before taking this seriously as a determination of 
the A parameter, we have to consider systematic errors. 

The most important question is whether the A-parameter extracted from a two-loop 
fit with the bare coupling is stable when going to three-loops, and when considering 
other schemes for the coupling. It is believed that perturbative scaling holds better 
when using an "improved" coupling The precision data for the string tension that 
are now available for a wide range of lattice spacings provide an excellent testing ground 
of these ideas. 



We will consider the E'-scheme [25], the ^/-scheme [26. ^J], and the MS scheme. In the 
-E-scheme (E for energy) one uses the leading perturbative expansion of the expectation 
value of the plaquette, Wn, to define g\ for SU(AT) gauge theory via 

2 4 2 N 2 -1 

\-W\\=w\g + w 2 g + ... = wig E , wi = — . (5.1) 

In practice one uses the value of the plaquette measured in a Monte Carlo simulation on 
the left hand side. 

To obtain the coupling in the V-scheme, originally defined via the static potential, 
we use the expansion of the logarithm of the plaquette to second order to first define 
the "plaquette coupling" g 2 P {q) |p4fl . Using the measured value of the plaquette, g 2 P {q) is 
obtained by solving, for pure SU(iV), 



lnWii = wigp(q) 



, . / 6.711706 \ , 
26 ln( ^ J g 2 P {q) 



(5.2) 
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Figure 4: Plot of the parameterization ( |4,4j ) of the string tension as a function of 
(3 = 6/<? 2 . Also shown are the data from table |] (o) that lead to the parameterization, 
as well as earlier results (x) for coarse lattices from refs. [||, 0- For clarity a few data 
points have been slightly displaced around their actual value. 
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Scheme 


cs 




bare 


0.189604 


28.80934 


E 


0.011629 


13.88031 


V 


-0.136923 


0.625192 


MS 


-0.012631 


1 



Table 5: Three-loop coefficient cs = (6f — ^o^2,s)/(2^q) of the scaling function and A- 
parameter for pure SU(3) gauge theory in various schemes. 



The optimal scale q to extract the coupling, according to [24], is given by the value 
q* = 3.4018/a M, M. 



By construction gp = g"y + O(gy). The exact coefficient of the third order term 
is now also known (cf. [^7], |2^, |2£| with the last missing ingredient recently provided 
in so that we can extract g v to third order from gp defined above:[] 



i + 



2.814089 2 t 
(4vr) 2 



■9W 



(5.3) 



for pure SU(3) and with q* as above. We can also obtain g^ to third order from gp, 



for which we will use the BLM [32] improved relation 

,2 / -5/6„*\ 2 / „* 



9 Z p(<1* 



1 2 . *.o 0.9546933 2 . *. 4 

i + ^ 9 W? + (47r)2 al(<f) A 



(5.4) 



In all four schemes the three-loop running of the coupling is known (see [| 
for the bare, e.g. [^| for the E-, [^] for the V-, and [31] for the MS-scheme). The 
scheme-dependent three-loop correction to the two-loop scaling function amounts to a 
multiplicative factor of the form 1 + csg^- The value of the coefficient cs in various 
schemes is given in table ||, together with the A-parameters. Our ansatz for a three-loop 
fit of the string tension is therefore of the form 

(aV5)M?)) = (1 + csgUq))f(g S (q)) [l + c 2 a{g s {q)f + c 4 a{g s {q)f] . 

(5-5) 

Ultimately, once better data for ro are available on fine lattices, these fits should also be 
performed for ro- At present, however, the data for the string tension are significantly 
better, overall. 

For each scheme we have performed fits using a two- loop and a three- loop ansatz. 
For the V and MS cases the coupling is then also extracted consistently to second, 
respectively, third order from the measured value of the plaquette. We find that all data 

7 Since the third order expansion of the plaquette is now also known in terms of gy, one could try to 
obtain the latter directly from the plaquette, without the detour via gp. However, it turns out that the 
cubic equation for g v has no physical, i.e. real and positive, solutions for (3 < 5.7. The detour via gp is 
therefore to be preferred. 
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bare 
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6/0 


l 


0.171(17) 


0.01325(12) 


178(2) 




3 




i 


0.157(16) 


0.01550(14) 


208(2j 


E 


2 


3(1 - Wu) 


l 


0.080(6) 


0.03817(21) 


246(1) 




3 


3(1 - W u ) 


l 


0.078(6) 


0.03864(21) 


249(1) 


V 


2 


9 2 p(q) 


3.4018 


0.057(6) 


0.2604(13) 


258(1) 




3 


eq. Q 


3.4018 


0.130(8) 


0.2449(15) 


242(1) 


MS 


2 


^(9*) + 3F^(ff*) 


1.4784 


0.074(6) 


0.3539(19) 


243(1) 




3 


eq. Q 


1.4784 


0.072(6) 


0.3647(19) 


251(1) 



Table 6: Fit results for the string tension data (3 > 5.6 in various schemes to an ansatz 
of the form (|5.5| ) (the values of the fit parameter C4 are not shown). The order indicates 
which scaling function was used (two- or three-loop) and, for the V- and MS-schemes, 
to which order the coupling was extracted from the measured plaquette. We assume 
Y / o : =465 MeV in converting to physical units. Details are in the text. 



can be fitted reasonably well in all schemes with either the two- or the three-loop ansatz. 
In either order the fits in the different schemes are of comparable quality. 

When leaving out small f3 points the fit parameters are essentially stable within errors 
in all cases. This is particularly impressive for As/iaq^/a), which is stable on the 1% 
level even when all (3 < 6.0 points are left out of the fit. By then C4 has become ill- 
determined and can be set to 0. But C2 is still very significantly different from 0; fixing it 
to gives very small confidence levels (often astronomically small; in the bare coupling 
scheme this is even true when leaving out the /3 = 6.0 points!). 

The stability of the fit parameters also holds when leaving out large j3 points. The 
A-parameter is stable on the 1% level even when leaving out all (3 = 6.2, 6.4 and 6.5 
points! This clearly illustrates the quality of the fit ansatz. 

We should mention that the reduced x 2 of the fits become very small (x 2 /^df < 0.25) 
once only (3 > 5.85 points are included in the fits. We take this as an indication that 
the errors on fine lattices have been, if anything, overestimated (at least for the data 
collected in table |j). 

In table |6| we show results for the various fits. For definiteness we have chosen to 
fit (3 = 5.6 — 6.5 in all cases. We now see interesting differences between the schemes: 
Most obviously, that the bare scheme has the largest a 2 coefficient C2, and shows a large 
increase in the A-parameter when going from two to three loops.0 The E- and MS- 
schemes, on the other hand, show much smaller a 2 corrections, and the A-parameter is 
much more stable when going from two to three loops. We regard the results in table ^ 

8 Note that the a 2 corrections to the perturbative scaling of the string tension are not "physical", 
in the way that, say, 0(a 2 ) violations of rotational symmetry would be. In fact, we could perform a 
(non-perturbative) change of coupling to the "cr-scheme", by defining g 2 . via (a^/a)(g) = f (g^) y/a / 'A. 
In this scheme there would be no 0(a 2 ) corrections whatsoever. 
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as a reconciliation of the ideas of Parisi-Lepage-Mackenzie with those of Allton. 

Note that the change in going from two to three loops is significantly larger for the 
V- than for the E- or MS-schemes. This can be traced back to the large value of the 
three-loop coefficient in the scaling function, cf. table [|. This is also the reason for the 
jump in C2- What is a bit surprising, perhaps, is that the results labeled MS, which 
really come from a mix of the MS- and the U-scheme (the latter was used to extract a 
physical coupling from the measured plaquette), show much better agreement between 
the two and three loops results than the U-scheme. In any case, one thing does seem 
clear, namely, that the values of A^jg extracted from the bare coupling scheme should 
not be trusted. 

As another check of systematic errors we have considered using the 1x2 Wilson 



loop, W\2i instead of Wu, to extract gp(q*)- Using the second order results of 33] 
extrapolated to infinite volume and q* = 3.0668/a in a two-loop fit, we find that A^jg 
goes up by 16 MeV and 14 MeV, respectively, compared to the results of table |6] for 
the V- and MS schemes. One can also define an "£^i2"-scheme, using the expansion 



of 1 — W\2 analogously to eq. (^TJ). From Q we obtain Aj^g/A^ 12 =27.4576 for this 
scheme. Two-loop fits now show a decrease of A^jg by 15 MeV compared to table ^. 
Unfortunately we can not perform three-loop checks for the "12" schemes. 

Our final check of systematic errors consists of the following worst case scenario. Let 
us assume that for some unknown reason the true ^Ja at both /3 = 6.4 and 6.5 are higher 
or lower by 2% than the values in table First of all we find that the x 2 °f our fits jump 
up by a factor of 2 or 3 (depending on the fitting interval), making the scenario appear 
unlikely indeed. Nevertheless, the value of Aj^g changes only by about 1 — 2%. This was 
to be expected, since, as mentioned earlier, we can leave out the (3 = 6.4, 6.5 (and even 
6.2) points completely from the fit, without significantly affecting the fit parameters. 

As our central value for A^jg we will quote the average of the three-loop results in the 
improved coupling schemes. Among these schemes, the systematic errors appear to be 
largest in the U-scheme; the difference between the two- and three-loop results, as well 
as switching from Wu to W±2 to extract the coupling, give a change in Aj^g of 16 MeV. 
We will quote this as our error: 

A m = 247(16) MeV . (5.6) 

This is consistent with the value recently obtained by the ALPHA collaboration |34j , as 
well as earlier results obtained from the string tension (see e.g. |l6|]). 

To get ([5j]) we used a value of \fo = 465 MeV for the string tension. This value 
is somewhat higher than the ones commonly used in the past. However, several recent 
studies quote a value of this size |(| , and we would also like to make the following 
simple, numerical observation. The results in table || show that r^y/a is consistent, 
within errors of about 1%, with the string picture value 



ro^/a = Jl.65 - ^ = 1.178219. (5.7) 

It is generally believed Q that tq is not larger than about 0.5 fm. Using ^fa = 1.178/Vo 
therefore gives 465 MeV (as a lower bound) for yfa. 
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6 Conclusions 



We have presented a high statistics study of the static potential in pure lattice gauge 
theory, using the Wilson plaquette action. These data and improved analysis strategies 
have enabled us to give very precise determinations of the string tensions and Sommer 
(-type) scales for (3 < 6.0, i.e. on lattices of spacings of about 0.1 fm and above. 
Combining our results with previous determinations, a\fa is now known to 1% in the 
range 5.54 < (3 < 6.5. According to our and other studies (quoted in sect. 4) this error 
also includes possible finite volume effects. 

The scale ro is presently known somewhat less accurately on the finest lattices. Using 
our method to determine Sommer-type scales it should be straightforward to determine 
ro for all lattice spacings at least as accurately as the string tension. In particular, 
any remaining finite volume effects should be much smaller for ro- Needless to say, our 
method can also be used for other kinds of gauge actions. 

We have employed the ansatz of Allton [17] to provide parameterizations of ^fa and 
ro that are accurate to about 1% for 5.6 < f3 < 6.5. By considering fits of this kind 
for the string tension in various coupling schemes we have estimated the A-parameter 
of quenched QCD to be A^jg = 247(16) MeV. We concluded that when it comes 
to extracting physical parameters, instead of just fitting data, the improved coupling 
schemes are very much superior to the bare coupling scheme, even when corrected for a 2 
effects as suggested by Allton. 

The results presented here should be useful in high precision scaling checks, 
comparisons of different actions, etc. Our motivation for this study has, in fact, been the 
desire to provide a significant check of the Symanzik improvement program by performing 
a scaling check of the rho mass in units of the string tension for the SW action with a non- 
perturbatively determined clover coefficient ||]. This is not possible without accurate 
knowledge of the string tension (or ro). Previous estimates of the string tension would 
have lead to misleading conclusions about the scaling behavior of the rho mass or other 
observables. 
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